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Abstract 

We discuss the notions of resurgence, formalizability, and 
formation of singularities in the context of partial dif- 
ferential equations. The results show that Ecalle's how 
analyzability theory extends naturally to PDEs. 



1 Introduction 

The study of nonlinear partial differential equations in the com- 
plex domain and especially formation of spontaneous singulari- 
ties of their solutions is not a well understood subject. The the- 
ory of Ecalle's anaiyzabie functions, originally developed (mainly) 
for functions of one variable, provides a set of tools which are 
well suited to address some of these issues, but the extension to 
several variables is not immediate. 

In the case of linear ODEs under suitable assumptions, there 
is a complete system of formal solutions as transseries, ^T] , and 
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these are generalized Borel (multi)summable to a fundamental 
system of actual solutions of the system; sufficiently powerful 
results of a similar nature have been shown in the nonlinear case 
as well [TT], [T3], [IS], [TE]. This is true for difference equations 
as well, El El- 
While transseries solutions can be considered and their Borel 
summation shown in the context of PDEs there are a number of 
difficulties specific to several variables. We first discuss a number 
of specific obstacles to extending the theory in a straightforward 
way, and then refer to what we expect to be a general approach to 
many problems and overview a number of recent results utilizing 
this approach. 

2 Difficulties of formalizability and an- 
alyzability in PDEs 

2.1 Insufficiency of formal representations 

For PDEs even the notion a general formal solution appears to 
elude a definition that is reasonably simple and useful. 
Example 1. The equation ft + f x — has the general solution 

f(t,x) = F(t-x) 

with F any different iable function; it is not clear to us that a 
worthwhile definition can be associated to the description "gen- 
eral formal expression in t — x" ; on the other hand, restricting 
ourselves to special, well defined, combinations in t — x would 
correspondingly limit the number of associated actual functions, 
precluding a complete solution of the original PDE. In this ex- 
ample, actual solutions outnumber by far formalizable ones. 

The (apparently) opposite situation is possible as well. 
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2.2 Absence of general summation procedures 

Let now /o G C°°(a,b). The initial value problem 

/t + i/* = 0; /| t=0 = /o(a;) (1) 
always has formal solution as £ j 0: 

r-tt&fu (2) 

fc=0 

but it has no actual solution, if /q is real-valued non-analytic (cf. 
the proof of next proposition). 

There is no nontrivial summation procedure of formal Taylor 
series over an interval nor a more restricted one that would asso- 
ciate actual solutions to (JTJ as the following proposition shows. 
(See also Remark [TJ). 

Proposition 1 Let S be a summation procedure defined on a 
differential algebra D s of formal series of the form 

oo 

9 = ^9k(x)t k 

where S is assumed to have the following (natural) properties: 

• S is linear. 

• S commutes with differentiation. 

• S(g) ~ g (x) as t 10. 

• If f G D s then Sf : (a, b) x (0, e) i— > C (where e is allowed 
to depend on f ). 

Assume if fo is real valued. Then f in f^) is in Ds iff / is 
convergent (in the usual sense). 
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Proof. Assume / G D s . Then, by the properties of S, the 
function / = Sf is a differentiable function, and it is a solution 
of (P) in a domain V = (a, b) x (0, e). If we write / = u + iv we 
see that the pair (u, v) satisfies the Cauchy-Riemann equations 
in T> and thus / is analytic in z = t + ix with (x, t) G T>. The 
third property of S shows that f is the limit as z approaches the 
interval (a, b) from the upper-half plane. Since f is real valued, 
then by the Schwarz reflection principle / extends analytically 
through (a, b) to a neighborhood of (a, b) in the lower half plane; 
in particular / is analytic on (a, b). But then (jTfl. which is the 
Taylor series of / at points on (a, b) is convergent. □ 

2.3 Obstacles to determining the formal solu- 



We now contrast formal analysis of ODEs and PDEs. 
1. Consider the Painleve PI equation 



a rather nontrivial example of a second order nonlinear 
ODE. A detailed analysis of transseries and their Borel 
summability in this example are discussed in ^Hj ■ We only 
mention a few aspects relevant to the present discussion. 

Finding formal solutions of is quite straightforward. 
Searching first for algebraic behavior, dominant balance 
shows that 6y 2 ~ —x, say y ~ 6~2txz, and then, consis- 
tent with this, y" = o{x). It follows that a formal series 
expansion can be gotten by taking y = and then, for 
neN, iterating the recurrence 



tions 



y = by + x 



(3) 



Vn+l 



6 2 2a/x — y: 
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A power series solution is readily obtained in this way, 



y = Q ^-^ + ^s^ + - (4) 

which is not classically convergent but is Borel summable 
to an actual solution [THj; the complete transseries can be 
calculated and Borel summed in a similar way, ^B]- The 
possibility (and convenience) of the formal calculation is 
partly due to asymptotic simplification, resulting in a dom- 
inant balance equation, 

6y 2 + x = 

which can be solved exactly, from which a complete solu- 
tion of the full problem follows by appropriate perturbation 
theory. 

2. Compare this problem to the periodically forced Schrodinger 
equation 

dib 

i = [-d xx + V(x) + Q{x) cos out] V (5) 

Under physically reasonable assumptions ip is transseriable 
® 

Theorem 2 ([3]) Assume Q,V are compactly supported 
and continuous, and Q > throughout the support of V. 
Fort > there exist N G N and {Tk}k<N, {Fuj ; k{t,x)}k<N, 
2tt / u -periodic functions oft, such that 

N 

i>(t,x) = z) + ^e- r *% fc (t,x) (6) 

j&L k=l 
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with 3?r fc > for all k< N, andhj{t,x) = O (H, 
have Borel summable power series in t, 

hj (t, x) = cbJ2 h kj (x)r k/2 (7) 

fc>3 

The operator LB (Laplace-Borel) stands for generalized 
Borel summation Ii6} /. 

Insofar as a formal analysis would be concerned, it is to be 
noted that there is no small parameter in (JHJ) and largeness 
of t does not make any term negligibly small; a posteri- 
ori, knowledge of the transseries (JHJ) confirms this. In this 
sense, (J5J) admits no further simplification. There is, to the 
knowledge of the authors, no straightforward formal way 
based on (J3J) to determine whether if) — > 0, let alone its 
asymptotic expansion. 

3 Overcoming these difficulties: the 
approach of asymptotic regulariza- 
tion 

First note an implication of Ecalle's analyzability techniques [S]: 
a wide class of problems can be regularized by suitable Borel 
transforms. Summability of general solutions of ODEs or dif- 
ference equations, HH1 EH HE] shows that, under appropriate 
transformations, the resulting equations admit convergent solu- 
tions, an indication of the regularity of the associated equation. 

Transseries are obtained, by suitable inverse transforms, from 
these regularized solutions. 

In the case of PDEs it appears that regularizing the equation is 
in many cases the adequate approach. The result © is obtained 
in this way. 
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3.0.1 Elementary illustration: regularizing the heat equa- 
tion 

/** - ft = (8) 
Since flSJ) is parabolic, power series solutions 

fe=0 fc=0 

are divergent even if Fo is analytic (but not entire). Nevertheless, 
under suitable assumptions, Borel summability results of such 
formal solutions have been shown by Lutz, Miyake, and Schafke 
[TUj and more general results of multisummability of linear PDEs 
have been obtained by Balser j7j. 

• The heat equation can be regularized by a suitable Borel 
summation. The divergence implied, under analyticity assump- 
tions, by (jOJ is Fk = 0{k\) which indicates Borel summation 
with respect to t~ x . Indeed, the substitution 

t = l/r; f(t,x)=t-^ 2 g(r,x) (10) 

yields 

2 1 
9xx + r g x + -rg = 

which becomes after formal inverse Laplace transform (Borel 
transform) in r, 

3 

P9pp + -j9 P ~ 9xx = (11) 

which is brought, by the substitution g(p, x) = p~^u(x, 2p? ); y = 
2p2 , to the wave equation, which is hyperbolic, thus regular 

UXX - Uyy = 0. (12) 

Existence and uniqueness of solutions to regular equations is 
guaranteed by Cauchy-Kowalevsky theory. For this simple equa- 
tion the general solution is certainly available in explicit form: 
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u = — y) + u + (x + y) with arbitrary C 2 functions. 

Since the solution of (fT2|) is related to a solution of (jSJ) through 
(fTUjl . to ensure that we do get a solution it is easy to check that 
we need to choose u_ = u + = u (up to an irrelevant additive 
constant which can be absorbed into tt_) which yields, 



y 



f(t,x) = t 2 y y-2 Mo (x + 2^) +m - 2 e -»/*d ; 

(13) 

which, after splitting the integral and making the substitutions 
x ± 2 y a = s is transformed into the usual heat kernel solution, 

f(t, x) = H «o(a) exp (-^^) ^ (14) 

In conclusion although there is perhaps no systematic way to 
formalize the general solution of the heat equation, appropri- 
ate inverse Laplace transforms allow us a complete solution of 
the problem (in an appropriate class of initial conditions which 
ensure convergence of the integrals). 

Remark 1 Proposition^ can be also understood in the following 
way. Equation (Qp is already regular. Any actual solution, if 
it exists with the initial condition given in the Proposition, is 
trivially formalizable since it is then analytic. It is thus natural 
that no further summable formal solutions exist. 

3.1 Nonlinear equations: regularizat ion by In- 
verse Laplace Transform 

In this section we briefly mention a number of our results of that 
substantiate regularizability. 
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3.1.1 

Consider the third order scalar evolution PDE: 

3 

ft~ fy m = ^b j (y,t;f)f^ + r(y,t); f(y,0) = fi (15) 

j=0 

Formal Inverse Laplace Transform with respect to y gives 1 

F t + P 3 F= Yl l B i* * (P' F ) * F * k ] + R & *) ( 16 ) 

j<3,k<oo 

where convolution is the Laplace one, (/ * g)(p) = f f(s)g(p — 
s)ds. This equation is regular in that formal power series in p 
converges, since the coefficients in the equation are analytic. 

Multiplying by the integrating factor of the l.h.s. and inte- 
grating yields 

F{p,t)=MF{p,t) 
= F Q (p, t)+ j (-l) j e- p3 ^ [&F) * B jik * F* k ] (p, r)dr 

j<3;fc<oo J ° 

The regularity of this equation plays a crucial role in the proofs 
in P where we find the actual solutions of equation (fT5]). 

3.1.2 

Similar methods were later extended 3] to equations of the form 

u t + V(di)u + g (x, t, {din}) = 0; u(x, 0) = u 7 (x) 

with u G C r , for t G (0, T) and large |x| in a poly-sector S in C d 
(dj. = d x \d x \...d x d d and j\ + ...+jd < n). The principal part of the 

'■For technical convenience, in pQ we used oversummation. The paper P] 
shows that in fact Borel summability holds in the correct variable, in the 
more general setting decribed in the next section. 
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constant coefficient n-th order differential operator V is subject 
to a cone condition. The nonlinearity g and the functions uj 
and u satisfy analyticity and decay assumptions in S. 

The paper [3] shows existence and uniqueness of the solution 
of this problem and finds its asymptotic behavior for large |x|. 

Under further regularity conditions on g and uj which ensure 
the existence of a formal asymptotic series solution for large |x| 
to the problem, we prove its Borel summability to the actual 
solution u. 

In special cases motivated by applications we show how the 
method can be adapted to obtain short-time existence, unique- 
ness and asymptotic behavior for small t, of sectorially analytic 
solutions, without size restriction on the space variable. 

3.2 Nonlinear Stokes phenomena and mov- 
able singularities 

In the context of ODEs it was shown ^5], under fairly general 
assumptions, that the information contained in the regularized 
problem (equivalently, in the transseries) can be used to deter- 
mine more global behavior of solutions of nonlinear equations, in 
particular the fact that they form spontaneous singularity close 
to anti-Stokes lines. The method, transasymptotic matching, 
was extended to difference equations [HI 13 El • 

In nonlinear partial differential equations, formation of singu- 
larities is a very important phenomenon but no general methods 
to address this issue existed. 

The method of regularization that we described provides such 
a method. We briefly discuss the main points of [2|. 

At present our methods apply to nonlinear evolution PDEs 
with one space variable; even for these, substantial new difficul- 
ties arise with respect to [TH] . 

Consider the modified Harry Dym equation (arising in Hele- 
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Shaw dynamics) 



H t -H 3 H xxx + H x -}-H 3 = 0; H(x, 0) = 

2 Jx 



in an appropriate sector. 

Small time behavior. From pQ it follows that there exists a 
unique solution to above problem, and it has Borel summable 
series for small t and small y = x — t: 

-1/2 ( 15 1 \ 2 f 25875 195 3 , 

Singularity manifolds near anti-Stokes lines. To apply the 
method of transasymptotic matching, we look on a scale where 
the asymptotic expansion becomes formally invalid: y = x — t = 
0(t 2 / 9 ). The transition variable is thus 

V = ^, r = t"\ H(x,t)=t-^G( V ,r) 

Substituting into (JT7J), we obtain the following equivalent equa- 
tion 

-rjG v + -tG t + -G 3 - G 3 G nm = 

9 9 ' v 9 2 mv 

The natural formal expansion solution in this regime is 

oo 

G(r ] ,r) = J2r k Gk(v) (17) 

fc=0 

with matching conditions at large r], to ensure the solution agrees 
with the one obtained in pQ: 

G { V ) ~ v - 1/2 ; GM M 



^k+l/2 



We show that the series (J 17)) is actually convergent and equals 
H{x, t) in the Borel summed region (the radius of convergence 



11 



shrinks however with if). The convergence problem is subtle 
and required a rather delicate construction of suitable invariant 
domains. Having shown that, it is intuitively clear (and not 
difficult to prove) that if Go is singular, then H is singular. The 
leading order solution Go satisfies 

G + 2 7 ] G' + 9G 3 G' " = 

while for k > 1, 

G CkGk = Rk 

where 

C k u = u + —r,u-^— + —j u 
and the right hand side Rk is given by 
Rh(v) = 2 G kl G k2 Gk 3 + ^2 G kl Gk 2 G k3 G'l! 4 

k 1 +k 2 +k 3 =k-l kj<k,J2kj=k 



The nonlinear ODE of Go with asymptotic condition has been 
studied in [T7j and computational evidence suggested clusters of 

1/3 

singularities rj s , where G (r]) ~ e 7 ™/ 3 (y) (77 — r] s ) 2 ^ 3 - 

G + 2 7/G , + 9G i G " = (18) 

For a rigorous singularity analysis of (fTSj) we now used transasymp- 
totic matching as developed for ODEs ^S]- The asymptotic be- 
havior of Go is of the form 

G (v) ~ V- 1/2 U(C) 

where 

q ■id\/2 
C = -\nC+ -In 77 + -i-rf' A + (2n - l)m 
8 27 

(where G is the Stokes constant) and U(() satisfies algebraic 
equation 



l e C + 2 = e -2VU (VU + 1 

Vu-i 
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Singularities of U(() occur at ( s = — In 4+ 2 — Z7r, corresponding 
for neN large to 

?4\/2 , Q 

— — rfJ A + - In = -2 + In 4 - (2n - l)nr + InC 
27 8 

The Theorem that we prove in [2j is that For a singularity fj s 
of G 2 there exists a domain V around the singularity 3 such 
that the expansion is convergent for small r. 

In particular, for small r, the singularity of G(rj, r) = H(x, t) 
approaches the singularity of Go and is, to the leading order, of 
the same type, (77 — r] s ) 2 ^. 
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